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I.  INTRODUCTION 


In  recent  years  there  has  been  a  great  interest  in  the  development 
of  sources  of  high  power  millimeter  waves  for  both  electron  cyclotron 
resonance  heating  of  fusion  plasmas  and  millimeter  radar.  Perhaps  the 
most  promising  of  these  sources  is  the  gyrotron,  where  the  rotational 
energy  of  a  spinning  electron  beam  is  coupled  into  a  RF  wave.  The  rela¬ 
tive  simplicity  of  fabrication  and  high  theoretical  efficiency  (>50  per¬ 
cent)  have  made  gyrotrons  particularly  attractive  for  these  applications. 

However,  regardless  of  these  advantages  the  gyrotron  is  nonetheless 
a  vacuum  tube  device  and  as  such  requires  an  output  coupling  window  in 
order  to  allow  extraction  of  energy  while  maintaining  vacuum  integrity. 
Should  this  window  present  a  poor  match,  a  significant  portion  of  the 
output  power  of  the  tube  will,  through  reflection  and/or  loss,  be  pre¬ 
vented  from  being  transmitted,  decreasing  the  actual  efficiency  of  the 
tube . 


Current  methods  of  window  design  involve  the  use  of  lumped  circuit 
element  equivalents  to  model  the  window  region.  These  techniques  assume 
that  the  region  in  which  the  window  is  located  is  cutoff  to  all  but  one 
(or  at  most  very  few)  modes.  The  physical  dimensions  of  the  window  are 
then  determined  pursuant  to  the  desired  RF  characteristics. 


These  methods  are  not  overly  accurate  when  applied  to  windows  to  be 
used  in  gyrotrons.  Unlike  most  other  microwave  tubes,  a  gyrotron  produces 


*  high  frequency,  high  power  RF  lignel  propagating  in  a  highly  overnoded 
waveguide.  Therefore,  a  technique  for  designing  gyrotron  windows  is 
required. 

This  report  presents  one  way  in  which  such  a  design  could  be  accom¬ 
plished.  It  is  based  on  methods  developed  in  the  field  of  laser  optics, 
a  regime  noted  for  severely  overmoded  high  power  systems.  In  addition, 
this  paper  contains  a  routine  which  allows  for  the  estimation  of  the 
mode  content  of  the  energy  intercepting  upon  the  window. 


II.  THEORETICAL  FORMULATION 


A.  MODE  CONVERSION  IN  TAPERED  WAVEGUIDES 

In  order  Co  accurately  design  a  window  for  maximum  efficiency,  it 
is  necessary  Co  ensure  that  Che  window  is  able  Co  transmit  all  of  Che 
modes  impinging  upon  it.  Whereas  a  gyrotron  cavity  is  designed  to  couple 
energy  into  one  selected  mode,  the  remainder  of  the  tube  geometry  is 
under  no  such  constraint.  As  can  be  seen  in  Figure  1,  there  is  an 
uptaper  from  the  cavity  into  the  collector,  and  a  further  dovntaper  from 
the  collector  to  the  window. 


The  net  effect  of  these  tapers  is  to  allow  energy  to  couple  out  of 
the  main  mode  into  undesired  spurious  modes.  Due  to  the  increased  diam¬ 
eter  of  the  tube  at  this  point,  these  spurious  modes  are  not  cutoff. 

Such  a  region  can  in  fact  allow  a  large  number  of  modes  to  exist.  As  an 
example,  at  60  GHz  a  region  2  inches  in  diameter  would  permit  263  modes 
to  propagate. 


The  magnitude  of  the  mode  conversion  can  be  computed  by  the  use  of 
the  generalized  telegraphist's  equations: 
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The  subscripts  i  and  p  indicate  various  modes.  Solymar  has  defined  the 
amplitude  of  the  forward  and  backward  waves  (A*  and  )  as 
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Rewriting  (1)  and  (2)  in  terms  of  aT  and  A.,  one  obtains 
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The  coupling  coefficient  from  the  pC  to  the  iC  mode  S1^  (  +  indi¬ 


cates  forward  coupling,  -  backward)  may  be  defined  in  terms  of  the  trans 


fer  function  T  .  of  (1,  2)  as 
Pi 


(7 


Solymar  has  provided  a  table  of  the  various  values  of  the  coupling  coef¬ 
ficient.  The  portion  of  the  table  of  interest  to  this  paper  is  repro¬ 
duced  in  Figure  2. 


The  majority  of  gyrotrons  in  production  today  operate  in  a  TE  mode 
so  equations  (5)  and  (6)  of  Figure  2  need  to  be  examined  with  greater 
scrutiny.  First,  define  the  axisymmetric  elemetary  wave  function  (the 
solution  to  the  Helmholtz  equation  in  a  cylindrical  axisynmetric  condue 
tor)  for  TE  modes  as 
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where  P  .  is  the  i  root  of  J'(x)  *  0.  However,  in  order  to  determine 
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can  be  seen  from  Figure  3,  the  normal  vector  n  is  identical  to  -  for  a 

cylindrical  waveguide.  Therefore,  one  finds  that 
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Applying  the  Bessel  function  derivative  relationship 
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Differentiating  again  with  respect  to  c  yields 


2  .  — n£-_  .  J  j  cos  n  ;  e-jk*z 


Let  us  now  define  the  above  mode  Co  be  t  .  rather  than  w,  and  fur- 
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thermore  let  there  exist  some  other  wave  function  i  .  such  that 
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Upon  examination  of  the  formula  for  the  coupling  coefficient,  it  can  be 
seen  that  the  critical  terms  consist  of  two  similar  line  integrals  around 
the  waveguide.  Due  to  this  similarity,  it  will  be  necessary  to  examine 
only  one  of  these  integrals,  allowing  the  complimentary  solution  to  be 
written  from  duality. 


The  integral  to  be  evaluated  is 
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This  equation  can  be  simplified  to 
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which  is  easily  solved  through  application  of  the  well  known  orthogonality 


condition  for  trigonometric  functions 
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Equation  (13)  can  therefore  be  expressed  at 
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A  similar  expression,  differing  only  in  subscripts,  is  obtained  for  the 
other  line  integral. 


The  necessary  calculations  having  been  performed,  it  is  now  possible 

to  write  the  closed  fora  solution  for  the  coupling  coefficient  from  one 

TE  mode  to  another.  Letting  represent  the  original  (pC^)  mode,  and 

.  the  mode  into  vhich  power  is  coupled,  one  obtains 
BJ 
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The  same  methodology  can  be  used  to  derive  the  coupling  coefficient 
from  TE  to  TM  modes.  Defining  the  TM  wave  function  to  be 
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it  can  be  shown  that 
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When  equations  (18)  and  (19)  are  introduced  into  the  formula  for 
the  coupling  coefficient,  one  obtains  an  equation  of  the  form 
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ip 


f 


cos  m$  sin  n$  dp 


There  remains  only  to  define  an  expression  for  the  backward  self 
coupling  term  S7^.  Again,  following  Solymar's  equations,  this  term 
involves  the  evaluation  of 
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Equation  (19)  has  provided  an  expression  for  3*/3s,  and  noting  that 
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S7.  to  be  written  as 
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By  drawing  upon  the  above  relationships  (5)  and  (6),  it  is  possible 
to  obtain  the  amplitude  of  the  spurious  modes  caused  by  changes  in  geom¬ 
etry  of  the  tube.  It  is  possible  at  this  point  to  make  a  few  generalized 
observations  about  this  sef  of  modified  telegraphist's  equations: 


1.  There  will  be  no  coupling  between  modes  in  a  waveguide  with 
constant  radius  (tan  6  ■  dlnK/dz  *  0) 


2.  If  the  waveguide  is  axisymmetric,  no  energy  is  coupled  between 
TE  and  TM  modes.  However,  should  the  conductor  be  physically 
perturbed,  coupling  will  occur. 


3.  These  equations  are  valid  for  any  slowly  tapering  waveguide. 
For  this  paper,  only  a  linear  taper  was  considered,  but  more 
complex  geometries  can  be  handled. 


4.  In  order  for  coupling  between  TE  modes  to  occur,  the  azimuthal 
mode  numbers  of  the  two  modes  must  be  the  same. 


A  schematic  representation  of  mode  conversion  is  provided  in  Figure  4. 


B.  WINDOW  CHARACTERISTICS 


Before  an  optimization  procedure  can  be  devised,  it  is  first  neces¬ 
sary  to  accurately  define  the  RF  characteristics  of  a  ceramic  window. 

In  order  to  arrive  at  such  a  definition,  the  window  geometry  must  be 
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considered.  Figure  5  shows  a  representation  of  a  non-undercut  five  disc 
window  of  constant  diameter.  Three  different  dielectrics  are  used  in 
this  window,  and  the  left  face  is  exposed  to  vacuum  while  the  right  is 
exposed  to  air. 

In  order  to  analyze  this  window  design,  it  was  decided  to  use  the 

2 

two  port  transmission  matrix  approach  as  developed  m  Ramo  .  In  this 
method,  it  is  required  that  each  of  the  planes  indicated  in  Figure  5  be 
treated  as  a  two  port  network.  This  multiplicity  of  planes,  upon  closer 
examination,  resolves  itself  into  two  basic  types  of  plane  -  an  interface 
between  dielectrics  and  a  length  of  connecting  waveguide  (irregardless 
of  whether  air  or  dielectric  filled). 

It  is  perhaps  easiest  to  begin  by  defining  the  transmission  (T) 
matrix  for  a  connecting  waveguide  section.  It  is  obvious  that  this  will 
consist  of  a  phase  shift  betwen  the  forward  and  backward  ports,  or 
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These  matrices  are  used  at  planes  b,  d,  f,  h,  and  j. 


The  remaining  planes  constitute  dielectric  interfaces,  and  here  Che 
transmission  matrix  must  be  obtained  through  the  application  of  boundary 


Figure  5  Window  geometry. 


conditions.  It  is  known  that  across  a  dielectric  boundary,  the  tangential 
components  of  the  electric  and  magnetic  fields  are  continuous.  Referring 
to  Figure  6,  this  can  be  expressed  in  terms  of  the  incident  and  reflected 
wave  amplitudes  as 
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Rearranging,  one  obtains 
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Or,  expressing  (26)  in  matrix  form 
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The  final  transmission  matrix  for  the  system  is  the  product  of  the  indi¬ 


vidual  matrices,  or 


A  convenient  parameter  for  expressing  the  efficiency  of  a  window  is 


its  voltage  standing  wave  ratio  (VSWR),  the  ratio  of  the  minimum  and 
maximian  values  of  the  standing  wave.  The  VSWR  is  defined  as 
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where  o  is  the  reflection  coefficient  of  the  systen.  Returning  again  to 
the  two  port  representation,  it  can  be  seen  that  o,  being  the  ratio  of 
the  reflected  to  the  incident  wave,  is  equal  to  bj/Sj  when  a,  »  0. 
Therefore 
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Introducing  (30)  into  (29),  one  obtains 
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Equation  (31)  is  the  form  of  the  figure  of  merit  that  will  be  used  in 
the  subsequent  computer  program  for  window  analysis.  It  should  be  noted 
that  the  transmission  matrix  elements  used  are  complex  and  frequency 
dependent,  leading  to  a  finite  bandwidth  for  low  VSWR  transmission. 


An  additional  parameter  of  interest  to  window  designers  is  that  of 
power  loss  into  the  dielectric.  The  loss  can  be  accurately  represented 


By  using  the  series  expansion  e 


x 
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This  expression  will  be  used  in  Part  III  to  help  define  the  constraints 
on  window  optimization. 


All  of  the  calculations  needed  to  define  the  problem  of  mode  conver¬ 
sion  and  window  analysis  have  been  performed.  It  is  now  possible  to 
proceed  with  the  description  of  the  numerical  solutions  and  optimization 
technique . 


III.  NUMERICAL  FORMULATION  AND  OPTIMIZATION 


A.  SOLUTION  OF  THE  MODE  CONVERSION  EQUATIONS 

In  examining  the  coupled  differential  equations  (5,6)  for  coupled 
mode  amplitude,  it  must  be  remembered  that  in  order  to  allow  a  numerical 
solution  to  be  found,  the  boundary  conditions  must  be  defined.  Solymar 
has  stated  that  the  boundaries  are  as  follows: 

A*(0)  -  A  aT(0)  ■  0  A"(L>  -  0  A*(L)  -  0  (35) 

d  o  i  m  l 

where  the  subscript  m  denotes  the  original  (driving)  mode  and  i  the  cou¬ 
pled  spurious  modes.  These  conditions  can  be  interpreted  as  meaning 
that  the  system  is  fed  by  only  one  mode  and  that  no  modes  are  being  intro¬ 
duced  back  into  the  output. 

Upon  first  glance,  it  would  appear  that  the  ordinary  differential 
equations  plus  boundary  conditions  combine  to  yield  a  straightforward 
coupled  initial  value  problem.  A  closer  inspection  reveals  that  whereas 
the  initial  conditions  for  the  forward  inodes  are  given  at  z  =  0,  those 
for  the  backward  modes  are  given  at  z  ■  L.  It  would  appear  then  that 
the  standard  methods  of  solution  for  initial  value  problems  (Runge-Kutta , 
predictor-corrector,  etc.)  are  inapplicable,  and  that  the  system  of  equa¬ 
tions  must  be  modified  to  allow  solution  through  the  use  of  a  two  point 
boundary  value  technique,  such  as  shooting.  One  would  like  to  avoid 
such  a  modification,  as  it  greatly  complicates  the  problem,  requiring 
solution  of  a  higher  order  system  and  increasing  computing  time. 


3 

However,  an  alternative  ia  available.  SSad  has  indicated  that  a 

4 

technique  developed  by  Denaan  ,  known  as  invariant  imbedding,  has  proved 
applicable  to  this  problem.  Since  the  full  force  of  the  imbedding  tech¬ 
nique  was  not  required  for  the  problem  at  hand,  a  full  discussion  of  the 
invariant  imbedding  method  will  not  be  undertaken  here;  those  interested 
should  consult  Denman. 

In  order  to  apply  the  invariant  technique,  one  must  be  able  to  con¬ 
sider  the  set  of  equations 

^  -^Z-  »  a(z)  u(z)  +  b(z)  v(z)  +  e(z)  (36a) 


-  v(?I  „  c(2)  u(z)  +  ^(z)  v(z)  +  f(2)  (36b) 


u(0)  -  u  v(l)  -  0  (36c) 

o 

Inspection  of  Solyman's  equations  (5,6)  reveals  that  they  are  of 
this  fora,  if  the  following  substitutions  are  made: 


u(z) 

1 

M?V 

N 

v(z) 

-  A~(z) 

a  (z) 

*  -J  6i  »  d(z) 

b(z) 

1  d  lnKi 

2  dz 

-  -c(z) 

e(z) 

-  K  +  S'  A*) 

p\  ip  p  ip  p) 

f(r) 

-  IK  A+  + 

pi  ip  p 

s;  A-\ 
ip  p ) 

(37) 
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With  this  change  of  variables,  the  solution  may  proceed.  As  equations 
(36a)  and  (36b)  are  similar,  a  solution  will  be  presented  in  detail  only 
for  (36a). 


To  begin  with,  it  is  required  that  the  variable  u(z)  be  replaced 
through  the  use  of  a  modified  Ricatti  transformation: 

u(z)  •  r^(r)  v(z)  ♦  z^iz)  (38) 

where  r^(z)  and  r^(z)  are  unknown  functions  whose  form  will  be  found 
presently.  Differentiating  (38)  with  respect  to  z  yields 


du 

dz 


_  ,dv  ,  ,  « 

rl  dz  1  2 


Substituting  (36a)  and  (36b)  into  (39) 


(39) 


au  ♦  bv  ♦  e  *  r^(-cu  -  dv  -  f)  ♦  r^v  ♦  r^  (60] 

(Each  term  in  the  above  equation  is  a  function  of  z;  indication  of  depen¬ 
dency  has  been  dropped  in  the  interest  of  clarity.) 


Making  use  of  (38)  in  (40) 


a(rlV  +  r2)  +  b  v  +  e 


^-c(r^v  +  r2)  -d  v  -  f  J 


+  r '  v  +•  r' 
rl  2 


Collecting  terms 


v|(a  ♦  d)r^  ♦  b  ♦  cr 


(a  ♦  cr^rj  ♦ 


fti  *  *  '  rd  ' 


Finally,  aeCting  the  temi  in  brackets  to  zero 


r|  ■  b  +  (a  +  d)  r^  +  cr^ 


(41a) 


(41b) 


(a  +  cr^)  r2  fr^  +  e 


Considering  the  point  z  *  0,  it  can  be  seen  that  (38)  becomes 


u(0)  -  rx(0)  v(0)  +  r2(0) 


Letting  r^(0)  *  0  (as  it  is  an  arbitrary  function  its  initial  value  may 
also  be  arbitrarily  chosen),  it  becomes  apparent  that  u(0)  “  r^O)  •  uq. 


Similarly,  if  one  lets  v(0)  »  q^zMz)  ♦  q0(z),  one  obtains 


qj  *  (d  ♦  cr j )  qx 


qi  •  (f  ♦  cr.)  q 


2'  u 


(42a) 


(42b) 


with  q^(0)  ■  1  and  q^O)  ■  0.  Therefore 


q 2 (^)  VU>  +  q2(*) 


v(0)  -  q2(z) 


q:(z) 


q,U)  q2U)  “  q2(z) 

qj(z)  q^2) 


-.V. 

•  --v-w. 


^JTVV'’ 


imm 

7~s\y^ 

►  9  0*0  •; 

.-.-v-V 

.■»y 

*  ‘  V  '/  - 


*•  \ 

■  ■  V 


■  v 
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Pausing  now  to  look  back  over  the  praeeeding  equations,  it  can  be 
aeen  that  the  original  system  of  equations  (36a, b)  has  been  replaced  by 
a  set  of  four  differential  equations,  (41a, b)  and  (42a, b),  which  along 
with  their  boundary  conditions  form  a  coupled  initial  value  problmn. 

The  four  quantities  (r^,  r^,  q^,  qj)  can  be  related  beck  to  the  original 
u  and  v  through  the  use  of  (38)  and  (44).  The  relationship  between  these 
newly  defined  variables  and  the  quantities  expressed  in  Solymar's  equa¬ 
tions  is  diagrammed  in  Figure  7. 


It  remains,  however,  to  obtain  an  explicit  manerical  solution  to 
this  set  of  differential  equations.  To  do  so,  it  was  decided  to  use  a 
multistep  technique,  such  as  a  predictor-corrector,  rather  than  a  single- 
step  method  like  Runge-Kutta.  This  choice  was  based  on  the  fact  that 
single-step  methods  use  Che  values  obtained  at  only  the  previous  step  in 
order  to  calculate  the  function  value  at  the  subsequent  step,  whereas 
multistep  techniques  make  use  of  the  data  obtained  at  several  previous 
points  to  calculate  the  next  step.  Thus  the  multistep  method  is  better 
equipped  to  track  a  quickly  varying  function,  as  it  would  be  less  likely 
to  incorrectly  predict  the  next  function  value. 


The  predictor-corrector  method  carries  a  further  advantage  in  that 
it  utilizes  both  implicit  and  explicit  forms  to  obtain  a  solution.  An 
implicit  form  is  one  where  the  value  of  a  dependent  variable  at  a  subse¬ 
quent  step  is  defined  in  terms  of  itself,  i.e.,  yn+^  “  f(xn+^,  yn+j). 
However,  since  one  is  usually  interested  in  determining  y  , ,  this  method 


cannot  be  directly  applied.  Initially,  therefore,  it  is  necessary  to 


1  dgnKi  2 
“2  dz  rl 


1  dgnKi 

2  dz 


»v, 


1  dgnKi  \ 

+  2  dz  fVf2 


*2(siiV 

p 

+  V  fs  +  A+ 
\  ip  P 
p 


+  s. 

•p 


+  s 

ip 


/  ..  A  1  dgnKi  \ 

K  +2—  Vqi 


A 

P 


♦  S 

<P 


.  +\  ^  1  dgnKi 

AP  r?— 


q 


J 


q,  (0)  *=1  q2  (0)  «  0 

r,  (0)  »  0  r„  (0)  »  A 

1  2  0 


7  Riccatti  transforms  of  Solymar's  equations 


employ  Che  explicit  form  of  Che  solution,  one  where  the  subsequent  step 


is  defined  by  the  previous  step,  or  yR+j  ■  f(xn+i»  yn)* 


The  Adams-Moulton  predictor-corrector  used  to  solve  the  problem  of 
interest  uses  an  explicit  method  as  the  predictor  and  an  implicit  method 
(using  yn+1  obtained  from  Che  predictor)  to  correct  the  approximated 
solution.  Since  Adams-Moulton  is  a  fourth  order  method,  it  calculates 
the  new  function  value  from  the  three  previous  values.  In  order  to  obtain 
these  initial  points,  Runge-Kutta  was  applied  for  the  first  three  steps, 
whereupon  the  Adams-Moulton  form  was  used: 


Predictor 


(0) 


y3  +  [55  f(x3,y3)  *  59  f(x2,y2)  +  37  f(xi»yi)  "  9  f(xo,yo)] 


(45a) 


Corrector 


y4J)  *  y3  +  14  [9  f(x4’y40))  *  19  f(x3»y3>  '  5(x2,y2)  «■  fCxj.yj)] 


(45b) 


This  method  was  applied  to  the  system  of  equations  shown  in 
Figure  7  in  order  to  obtain  rj ,  r2>  q^ ,  and  q2  at  each  interval  in  the 
taper  region.  Once  these  values  have  been  found,  it  is  a  simple  matter 
to  obtain  the  magnitudes  of  the  forward  and  backward  waves  for  each  mode 
throughout  the  taper. 


The  analysis  described  above  was  performed  on  both  the  uptaper  and 
downtaper  regions  shown  in  Figure  1.  Since  the  coupling  coefficients 
are  equal  to  zero  in  regions  where  the  waveguide  diameter  is  unchanging, 
there  was  no  need  to  perform  calculations  through  the  actual  collector 


region  itself.  The  final  amplitudes  from  the  uptaper  were  used  as  the 
starting  conditions  for  the  downtaper. 


As  vas  mentioned  before,  only  a  linear  taper  vas  considered  during 
the  problem  solution.  However,  the  only  modification  needed  to  be  per¬ 
formed  to  the  above  routine  in  order  to  permit  the  handling  of  more  com¬ 
plex  geometries  is  the  addition  of  a  short  "preprocessor"  program  to 
supply  the  radius  and  can  6  values  at  each  iterative  step.  Also,  since 
this  project  dealt  with  window  design,  only  the  mode  content  at  the  win¬ 
dow  was  considered  important.  However,  this  method  will  allow  for  the 
solution  for  the  mode  content  at  any  point  in  the  collector  region,  as 
well  as  the  calculation  of  the  voltage  travelling  wave  ratio  (VTWR)  as 
defined  by  Stone^: 


*J  ♦  I*. 

where  A  and  A  are  the  amplitudes  of  two  modes.  This  VTWR  is  •  measur- 
n  m 

able  quantity  representing  the  beat  pattern  resulting  from  the  superposi 
tion  of  two  traveling  waves.  This  physical  parameter  allows  for  the 
possible  experimental  verification  of  the  mode  content. 


By  using  a  predictor-corrector  on  the  Iicatti  transformation  of 
Solymar's  modified  telegraphist's  equations,  it  has  become  possible  to 
characterise  the  mode  content  of  the  W  energy  incident  upon  the  window, 
and  can  now  proceed  to  the  actual  opt  imitation  of  the  window  thickness. 


B.  OPTIMIZATION  OF  WINDOW  THICKNESS 


When  designing  e  window,  it  is  desirable  to  have  the  lowest  possible 
VSWR  over  the  widest  possible  bandwidth.  One  implementation  of  this 
transmission  profile  is  the  Butterworth  or  maximally  flat  response,  where 
the  VSWR  increases  smoothly  from  unity  in  either  direction  from  the  center 
frequency.  This  response  function  is  desirable  in  that  minor  variations 
in  frequency  will-  not  cause  an  abrupt  change  in  transmission,  as  can 
occur  with  other  filter  (window)  responses. 

Achieving  such  a  transmission  characteristic  has  until  now  implied 
searching  for  the  combination  of  window  thicknesses  (assuning  known  diam¬ 
eters  and  dielectric  constants)  by  either  the  laborious  process  of  cold 
testing  various  combinations  or  repeating  the  calculation  of  the  window 
characteristic  with  varying  widths;  at  best,  a  tedious  procedure.  Whereas 
at  lower  frequencies  a  lumped  element  analysis  or  quarter  wavelength 
method  can  be  used,  these  both  assuse  that  the  impedance  of  the  window 
can  be  varied  to  exactly  match  the  calculated  value.  At  gyrotron  fre¬ 
quencies,  only  a  few  window  materials  are  available,  and  for  ease  of 
construction  it  is  desirable  to  obtain  a  constant  window  diameter,  so 
the  low  frequency  techniques  are  inapplicable. 

It  is  therefore  necessary  to  determine  another  method  of  optimiza¬ 
tion.  Optimization  by  itself  implies  a  search  for  the  best  condition, 
so  in  effect  the  problem  involves  seeking  a  combination  of  window  thick¬ 
nesses.  To  do  so,  a  technique  known  as  the  simplex  method  was  applied. 


The  simplex  Method,  as  developed  by  Nelder  and  Mead^,  is  a  means  of 
searching  for  either  maxima  or  minima.  A  simplex  is  an  area  in  n-space 
defined  by  n+1  vertices  which  are  linearly  indepedent  and  all  line  seg¬ 
ments  connecting  them.  The  vertices  consist  of  the  function  values  eval¬ 
uated  at  n+1  randomly  defined  points  in  n-space  (Figure  8).  This  rather 
odd  appearing  construct  locates  a  minimum  in  the  following  fashion: 
first,  the  point  associated  with  the  highest  function  value  P^  is  identi¬ 
fied.  This  point  is  replaced  by  a  new  point  defined  as  the  reflection 

of  fu  through  the  centroid  P  (the  center  of  all  points  other  than  f  ) 
n  n 

according  to 

P*  *  (1  ♦  a)  p  -  aP.  a  <  1  (47) 

n 

The  function  value  at  the  new  point  P*  is  now  found.  Should  it  prove  to 
be  the  lowest  function  value,  a  further  expansion  in  that  direction  is 
performed  by 

P**  -  YP*  ♦  (1  -  Y)  P  Y  >  1  (48) 

If  f(P**)  <  f. ,  then  F„  is  replaced  by  P**;  otherwise,  it  is  replaced  by 
L  n 

p*. 

Occasionally,  though,  it  will  oc-ur  that  the  reflected  point  P*  is 

associated  with  a  function  value  which  is  still  greater  than  that  at  any 

other  point,  i.e.,  f(P*)  ■  f  .  At  this  point  it  becomes  necessary  to 

n 

contract  the  point  min  (P„,  P*)  back  toward  the  centroid  using 


P**  -  8  ?„  ♦  (1  -  8)  P 

H 


o  <  e  <  l 


(49) 


If  the  problem  persists,  that  is  f(P**)  >  ,  it  is  an  indication  that  a 

minimizing  value  is  located  vithin  the  bounds  of  the  simplex,  so  the 
entire  simplex  (not  just  P  )  is  contracted  by  replacing  all  points  (P.) 

n  X 

by  1/2(P.  +  P  ).  This  shrinking  of  the  area  will  continue  until  the 

X  L 

minimum  is  found  or  is  placed  outside  of  the  simplex. 

A  simplex  can  be  envisioned  as  following  the  topology  steadily  down¬ 
ward,  its  speed  of  descent  in  direct  proportion  to  the  gradient;  the 
steeper  the  slope,  the  faster  the  simplex  proceeds.  Upon  reaching  a 
valley  containing  a  minimum,  it  contracts  about  it  until  the  stopping 
conditions  for  the  search  are  met.  The  simplex  method  does  however  have 
one  major  flaw.  It  is  unable  to  distinguish  between  a  local  and  a  global 
minimum  (Figure  9).  Once  it  enters  a  valley,  it  will  proceed  to  the 
lowest  point  in  the  valley,  regardless  of  whether  a  neighboring  valley 
contains  an  even  lower  function  value.  This  problem  will  be  addressed 
as  the  solution  develops. 

Before  beginning  the  search  itself,  it  proves  efficacious  to  define 
the  area  to  be  searched  a  bit  more  closely.  The  VSVR,  which  is  the  func¬ 
tion  to  be  minimized,  is  dependent  upon  the  thicknesses  of  the  n  windows 
and  the  frequency,  or  n+1  variables  altogether.  A  little  thought  reveals 
that  searching  all  n+1  variables  would  involve  looking  for  the  exact 
VSWR  vs  frequency  response  desired.  This  is  a  favorable  goal,  but  hardly 
a  practical  one.  Should  the  actual  available  response  differ  even 
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•lightly  from  the  ideal  values,  the  search  will  fail.  Additionally,  the 
•implex  technique  by  itself  performs  on  the  order  of  n  calculations  at 
each  step,  while  the  evaluation  of  the  function  (as  per  Section  II. A) 
employs  0(n  )  operations;  thus  it  would  be  beneficial  to  reduce  n. 

The  scheme  used  herein  avoids  both  these  problems.  Since  it  is 
known  that  the  lowest  VSWR  should  occur  at  the  center  of  the  band,  the 
search  used  in  this  project  held  its  frequency  variable  fixed  at  the 
center  frequency,  eliminating  that  variable.  Once  the  center  band  mini¬ 
mum  had  been  found,  the  response  over  the  bandwidth  would  be  examined. 

Up  until  now,  no  mention  had  been  made  of  exactly  what  area  is  to 
be  searched.  The  search  parameters  have  been  narrowed  down  to  only  the 
window  thicknesses,  but  searching  for  every  possible  width  combination 
would  prove  temporally  impossible  as  well  as  yielding  impractical  results 
(such  as  the  trivial  solution).  Two  limiting  factors  do  however  exist. 

The  first  is  associated  with  the  fact  that  the  window  must  bear  the  mechan¬ 
ical  stress  of  maintaining  a  vacuum,  so  that  it  must  be  thick  enough 
that  it  will  not  be  collapsed  by  the  pressure  difference.  At  the  same 

time,  it  will  be  remembered  that  the  power  absorbed  by  a  window  is 

1 

directly  related  to  its  thickness  (Eq.  34),  so  that  an  overly  wide  window 
would  absorb  a  large  mnount  of  RF  energy,  introducing  thermal  stress 
which  can  lead  to  cracking  and  window  failure.  These  two  mechanical 
considerations  therefore  define  the  upper  and  lower  parmneters  for  each 
of  the  search  variables.  In  practice,  they  can  be  determined  either  by 
a  detailed  mechanical  analysis  relating  maximum  permissible  power  loss 


and  stress  strength  to  thickness,  or,  ss  appears  to  be  more  commonly 
done,  by  intuitive  approximation. 


Nov  that  the  search  region  has  been  constrained,  one  further  obsta¬ 
cle  remains  to  be  removed.  It  vas  mentioned  earlier  that  the  simplex 
can  become  trapped  in  a  local  minimum  without  locating  a  global  minimum. 
This  problem  vas  circumvented  by  dividing  the  region  into  smaller  areas, 
thus  lowering  the  chances  that  one  area  would  contain  more  than  one  mini¬ 
mum.  A  simplex  vas  released  into  each  subdivision;  if  no  minimum  vas 
found  after  a  predetermined  number  of  iterations,  the  method  proceeded 
to  another  subdivision.  If  a  minimum  vas  found,  the  coordinates  of  the 
point  (window  lengths)  were  stored  and  searching  was  contained  in  another 
region.  The  stopping  criterion  used  was  that  the  function  value  should 
be  a  small  difference  away  from  the  minimum  value  (which,  for  VSWR,  is 
1 )  or 


VSWR  -  1  <  e  (50) 

After  every  subregion  has  been  examined,  each  discovered  minimum  is 
subjected  to  the  following  inspection  to  determine  if  it  meets  the  trans¬ 
mission  requirements.  The  value  of  the  VSWR  is  calculated  at  eight  addi¬ 
tional  points  -  center  frequency  (fc)  ♦  bandwidth/2  (BW/2),  fc  -BW/2, 
fc  ♦  BW/A,  fc  -  BW/A,  fc  ♦  BW/2  ♦  t,  fc  -  BW/2  ♦  t,  fc  ♦  BW/A  ♦  t,  fc  - 
BW/A  ♦  c ,  where  c  is  a  small  frequency  spacing  unrelated  to  the  c  in 
Eq.  50.  The  reason  for  calculating  the  function  value  at  pairs  of  closely 
spaced  points  is  to  enable  calculation  of  the  derivative  of  the  function 


at  these  points.  Zf  the  function  values  are  within  the  specified  liaits 
for  VSWR  across  the  band,  and  the  derivative  indicate  no  wild  oscillations, 
the  combination  of  thicknesses  is  declared  to  be  a  viable  one.  Occasion* 
ally,  aore  than  one  combination  will  be  passed,  whereupon  it  falls  to 
the  user  to  select  the  more  desirable  of  the  two,  whether  by  transmission, 
machining  ease,  or  lower  overall  power  absorbtion. 


IV.  CALCULATED  RESULTS 


It  is  always  expedient  (as  well  as  reassuring)  to  be  able  to  verify 
that  a  theory  coincides  with  actuality,  though  it  is  not  always  easy  to 
do  so.  It  was  not  within  the  scope  of  this  project  to  experimentally 
verify  the  accuracy  of  the  developed  design  technique,  so  in  order  to 
prove  its  validity  the  methods  proposed  in  this  paper  were  applied  to 
already  known  situations. 

To  confirm  the  mode  coupling  solution,  the  results  presented  by  the 
program  were  compared  against  those  obtained  by  a  similar  method  by 
Stone. ^  Having  access  to  greater  resources,  Stone  was  able  to  experi- 
mentally  determine  that  the  calculations  matched  reality  to  ±5Z.  There¬ 
fore,  agreement  with  Stone's  results  would  vindicate  the  method  proposed 
in  this  paper.  As  can  be  seen  in  Figure  10,  close  agreement  was  obtained. 
The  differences  that  occur  are  most  plausibly  explained  by  the  fact  that 
the  two  methods  used  different  routines  to  solve  the  differential 
equations . 

Window  optimization  was  validated  by  an  attempt  to  recreate  the 
values  currently  in  use  at  Hughes  Aircraft  Company.  The  three  frequen¬ 
cies  chosen  (28,  60  and  100  GHz)  represent  existing  gyrotrons  whose  win¬ 
dow  transmission  characteristics  are  considered  acceptable.  Shown  in 
Figures  11  to  16  are  the  parameters  used  in  the  calculations,  the  computer 
geneated  results,  and  the  actual  design  dimensions.  Discrepancies  in 
these  results,  though  slight,  ere  the  result  of  machining  considerations  - 


MODE  NUMBER 


-•  AFTER  STONE5 
■O  CALCULATED 


Figure  10  Mode  content 


DESIGN  PARAMETERS 


CENTER  FREQUENCY 

28  GHz 

PEAK  POWER 

100  kW 

BANDWIDTH 

4  GHz  BELOW  VSWR  -  1. 

MODE. 

TE  .  (p  .  ■  7.0155) 
o2  m 

NO.  OF  WINDOW  ELEMENTS 

3 

DIELECTRIC  CONSTANTS 

3.8.  9.6.  3.8 

LOSS  TANGENTS 

0.00158.0.00068.0.00158 

WINDOW  DIAMETER 

3  in. 

MAX.  POWER  LOST  IN 

EACH  WINDOW  ELEMENT 

0.5, 0.5, 0.5  kW 

MAXIMUM  SLOPE  (— V  "  ) 

0.35 

MAXIMUM  VSWR  AT 

CENTER  FREQUENCY 

1.05 

OPTIMIZED  WINDOW 

WINDOW  THICKNESS  (W) 

0.055. 0.0675. 0.055 

POWER  LOST  IN  EACH  SECTION 

0.13, 0.069, 0.13  kW 

Figure  11  28  GHz  window  parameters 


HUGHES  MULTI  DISC  WINDOW  DESIGN 


DISC  WIDTH  -0.0000 
WAVELENGTHS  •  0.0000 


CERAM.  WIDTH  -  0.0550  CEN  SPACING  -  0.0675 

WAVELENGTHS  -  0.2510  WAVELENGTH  -  0.4936 


DESIGN  PARAMETERS 


CENTER  FREQUENCY 

60  GHz 

PEAK  POWER 

200  kW 

BANOWIDTH 

3  GHz  BELOW  VSWR  - 

MODE 

TEo2  <Pni  *  7-01 56) 

NO.  OF  WINDOW  ELEMENTS 

5 

DIELECTRIC  CONSTANTS 

3.8,  9.6. 1.75.9.6,  3.8 

LOSS  TANGENTS 

0.00158.0.00068. 0.005. 
0.00068.  0.00158 

WINDOW  DIAMETER 

2  in. 

MAX.  POWER  LOST  IN 

EACH  WINDOW  ELEMENT 

1. 1.3. 1,  1  kW 

MAXIMUM  SLOPE  (— 

0.4 

MAXIMUM  VSWR  AT 

CENTER  FREQUENCY 

1.1 

OPTIMIZED  WINDOW 


WINDOW  THICKNESS  (IN)  0.073.  0.088.  0.0605. 

0  088.  0.073 


POWER  LOST  IN  EACH  SECTION  0.74, 0.39,  2.58,  0.39.  0.74  kW 


Figure  13  60  GHz  window  parameters 


HUGHES  MULTI  DISC  WINDOW  DESIGN 

DISC  WIDTH  -  0.0730  CERAM.  WIDTH  -  0.0880  CEN  SPACING  -  0.0805 

WAVELENGTHS -0.7188  WAVELENGTHS  -  1.3827  WAVELENGTHS  -  0.5339 


DESIGN  PARAMETERS 


CENTER  FREQUENCY 
PEAK  POWER 
BANDWIDTH 

MODE 

NO.  OF  WINDOW  ELEMENTS 

DIELECTRIC  CONSTANTS 
LOSS  TANGENTS 

WINDOW  DIAMETER 

MAX.  POWER  LOST  IN 
EACH  WINDOW  ELEMENT 

MAXIMUM  SLOPE  (-Vd^V-R) 

MAXIMUM  VSWR  AT 
CENTER  FREOUENCY 

OPTIMIZED  WINDOW 

WINDOW  THICKNESS  (IN.) 
POWER  LOST  IN  EACH  SECTION 


100  GHz 
500  kW 

5  GHz  BELOW  VSWR  -  1.5 
TE^  <Pn,  -  13.324) 

5 

3.8.86. 1.75. 9.6. 3.8 

0.00158, 0.00068. 0.005 
0.00068. 0.00158 

2  in. 

2. 2.  8.  2. 2  kw 

0.4 

1.05 


0.045.  0.052,  0.049. 0.052.  0.045 
1.90, 0.952, 6.52, 0.952. 1.90  kW 


Figure  15  100  GHz  window  perimeters 


HUGHES  MULTI  DISC  WINDOW  DESIGN 

DISC  WIDTH -0.0450  CERAM.  WIDTH  -  0.0520  CEN  SPACING  -  0.0490  . 

WAVELENGTHS -0.7371  WAVELENGTHS  -  1 .3607  WAVELENGTHS  -  0.5393 


FREQUENCY  (GHZ) 


Figure  16  100  GHz  window  response. 


it  is  impractical  to  specify  dimensions  to  tenths  of  mils.  An  additional 
note  -  it  was  experimentally  determined  that  slight  differences  in  window 
diameter  do  not  affect  the  VSWR  as  strongly  as  predicted  by  computer 
calculations. 

The  optimization  program,  being  written  as  part  of  a  research  pro¬ 
ject,  is  slightly  esoteric  in  its  operation  and  requires  a  slight  degree 
of  finesse  in  choosing  the  starting  conditions.  It  was  found  that  four 
factors  are  prisurily  responsible  for  those  numerous  occasions  during 
program  development  when  no  answers  were  returned.  First,  if  the  limits 
imposed  upon  the  lengths  are  too  close,  there  may  not  be  a  minimum  con¬ 
tained  in  the  region.  This  can  occur  if  one  attempts  to  specify  on  upper 
limit  (thermal  limit)  too  close  to  the  mechanical  limit,  i.e.,  constrain 
the  problem  severely  toward  the  thinnest  possible  window.  Another  prob¬ 
lem  arising  from  the  user  asking  for  too  much  comes  from  requiring  that 
the  VSWR  must  be  close  to  one  across  the  entire  band.  This  simply  cannot 
always  be  achieved. 

The  above  errors  are  caused  by  requiring  the  program  to  produce  an 
impossibly  ideal  design.  Two  further  factors  enter  as  a  result  of  the 
method  itself.  Should  the  subregions  be  made  too  large,  the  simplex  may 
trap  itself  in  a  local  minimum  and  be  unable  to  locate  an  optimum  solu¬ 
tion.  Similarly,  should  the  function  be  extremely  smooch,  the  simplex, 
whose  race  of  descent  is  determined  by  the  slope  of  the  function,  may 
not  be  able  to  converge  to  a  solution  within  the  maximum  number  of  stages 
These  two  probelms  can  only  be  avoided  by  developing  an  intuitive  feel 


V.  CONCLUSIONS 


As  was  shown  in  Che  proceeding  section ,  sgreenent  with  actual  results 
was  good.  Variations  occurred  either  from  different  calculation  tech¬ 
niques  or  alteration  of  results  to  meet  physically  realizable  conditions. 
This  similarity  of  results  leads  one  to  conclude  that  the  methods  presented 
in  this  paper  can  be  used  in  the  future  as  the  principle  method  of  window 
design. 

Originally,  it  was  intended  to  combine  the  mode  content  results 
with  the  window  optimization  routine,  designing  the  window  so  as  to  allow 
transmission  of  more  than  one  mode.  Unf ortunately,  however,  a  creditable 
means  of  weighing  the  modes  so  as  to  assign  priority  in  transmission 
characteristics  proved  unobtainable.  This,  along  with  allowing  for  curved 
window  surfaces,  are  fertile  ground  for  future  work. 


I 


BIBLIOGRAPHY 


1.  Solymar,  L.  Spurious  Mode  Generation  in  Nonuniform  Waveguide,  IRE 
Transactions  in  Microwave  Theory  and  T^hniquea,  July  1959,  pp.  379-383 

2.  Ramo,  S. ,  Vhinnery,  J. ,  Van  Duzer,  T.  Fields  and  Waves  in  Communica¬ 
tion  Electronics.  John  Wiley  1965. 

3.  Saad  S. ,  Davies,  J.  Computer  Analysis  of  Gradually  Tapered  Waveguide, 
IEEE  KTT,  May  1977,  pp.  437-440. 

4.  Deman,  E.  Coupled  Modes  in  Plasmas,  Elastic  Media,  and  Parametric 
Amplifiers.  Elsevier,  1970. 

5.  Stone,  D.  Measurement  of  VSWR  and  Mode  Content  in  Multimode  Wave¬ 
guides,  presented  Dec.  1,  1983,  Univ.  of  Utah  AFTER  Program  lecture. 

6.  Burden,  R.  Numerical  Analysis,  Prindle,  Weber  and  Schmidt,  1981. 

7.  Nelder,  J.  and  Mead,  R.  A  Simplex  Method  for  Function  Minimization. 


